A kinetic equation approach to the anomalous Hall effect in a diffusive Rashba 
two-dimensional electron system with magnetization 
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We present a two-band kinetic equation method to investigate the anomalous Hall effect in a 
Rashba two-dimensional electron system subjected to a homogeneous magnetization. The electron- 
impurity scattering is taken into account in the self-consistent Born approximation. It is demon- 
strated that the impurity-density-free anomalous Hall conductivity arises from an intrinsic and a 
disorder-mediated mechanisms, associated respectively with the electron states under and near the 
Fermi surface. The intrinsic mechanism relates to a dc-field-induced transition process, or in other 
words, a linear stationary Rabi process. The disorder-mediated one corresponds to a scattering 
between impurities and the electrons participating in longitudinal transport. Numerically, the de- 
pendencies of the anomalous Hall conductivity on the spin-orbit coupling constants and the strength 
of the magnetization are demonstrated for both the short- and long-range collisions. 

PACS numbers: 72.20.Fr, 73.50.Dn, 73.63.Hs 



I. INTRODUCTION 

In the presence of a magnetic field, the Hall effect in- 
duced by Lorentz force is a powerful tool for measurement 
of the concentration and the nature of the free carriers. 
However, in many ferromagnets, a nonvanishing trans- 
verse resistivity can also be produced by the spontaneous 
magnetization, which exists in the absence of external 
fieldsi This so-called anomalous Hall effect (AHE) has 
been extensively studied in ferromagnetic semiconductors 
in recent years £ 

AHE has been first predicted about five decades agoi^ 
At early stages, to interpret it, two mechanisms, namely 
skew scattering 4 and side-jump process^ have been 
proposed^ These two mechanisms are based on a spin- 
orbit coupling included into the potential of electron- 
impurity scattering. Obviously, the corresponding con- 
tributions to AHE should rely on the electron-impurity 
collision. 

The recent theoretical interest has been 
cused on another "intrinsic" mechanism 
AHEi 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 : 15 ! 16 ! 17 This mechanism 



been first discussed by Karplus and Luttinger,' 



fo- 
of 
has 
and 

clearly rederived by MacDonald et alm*^ and Nagaosa 
et alM^ By inclusion of the spin-orbit interaction into 
Hamiltonian of the free electrons, nontrivial transverse 
conductivity has been obtained even in the absence 
of disorder scattering. It makes clear that this AHE 
is associated with the Berry phase in the momentum 
spaced 

More recently, according to the well-known result of 
Stfeda in the context of Hall effect in magnetically two- 
dimensional (2D) systems^ Dugaev et ai. discussed an- 
other mechanism of AHE by means of Kubo formula. 20 
This mechanism is related to electron states in the 
vicinity of Fermi surface, and depends on the electron- 
impurity scattering but is independent of the impurity 
density iVj. By considering the short-range disorder col- 
lision, it has been found that this collision-related mecha- 



nism makes a contribution to anomalous Hall conductiv- 
ity of the same order of magnitude as the intrinsic one. 
However, as shown in Ref. l2fl the Kubo approach be- 
comes questionable when used to investigate the collision- 
related AHE. Taking the static limit uj — > before or 
after JV, — > leads to completely different results. At the 
same time, up to now, this AHE has been studied only 
for short-range electron-impurity scattering. We know 
that in realistic 2D semiconductors, the electron density 
is not large enough to screen the charged impurities. The 
Coulomb interaction between electron and impurity is in- 
evitably long-ranged. 

In this paper, we employ a two-band kinetic equa- 
tion approach to study the impurity-density-independent 
AHE in 2D Rashba electron systems with magnetiza- 
tion. Within this method, we can consider the time- 
independent nature of dc transport from the beginning 
and hence the problem caused by taking the static limit 
in Kubo formalism can be avoided. At the same time, 
the long-range electron-impurity collision can be easily 
handled. The obtained kinetic equations in the helicity 
basis allow us to interpret the above two mechanisms of 
AHE in terms of interband processes. The intrinsic AHE 
arises from a dc-field-induced direct transition between 
two unperturbed spin-orbit-coupled bands. In another 
point of view, this AHE can also be understood as a 
result of a quantum interference between perturbed elec- 
trons in different bands in the first order of dc field, i.e. 
as a linear stationary Rabi process. On the other hand, 
the electrons joining in the longitudinal transport can be 
resonantly scattered by the impurities, leading to another 
interband polarization independent of impurity density. 
This process corresponds to the collision-related mecha- 
nism of AHE. Similar picture has been demonstrated in 
the previous studies on the spin-Hall effect in 2D sys- 
tems with spin-orbit coupling^! Based on the derived 
equations, we numerically investigate the dependencies 
of the AHE on the spin-orbit coupling constant and the 
strength of exchange field. At the same time, the compar- 
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ison between the effects of short- and long-range disorders 
on AHE is made. 

The paper is organized as follows. In Sec. II we derive 
the kinetic equation for nonequilibrium distribution func- 
tions and attempt to resolve it. The above-mentioned 
two mechanisms of AHE correspond to two distinct com- 
ponents of the solution. In Sec. Ill we perform a nu- 
merical calculation to investigate the effect of spin-orbit 
coupling constant and magnetization on AHE. Finally, 
we conclude our results in Sec. IV. 



II. FORMALISM 

A. Kinetic equation 

We consider a quasi-two-dimensional system in the x — 
y plane with Rashba spin-orbit interaction. This system 
is subjected to a homogeneous magnetization Mo along 
the z-direction. The noninteracting Hamiltonian of the 
considered system has the form 



H Q = e p + a(a x p y - a y p x ) - Ma z , 



(1) 



where, p = (p x ,p y ) = (pcos(/> p ,psin(/>p) is the electron 
momentum, e p — p 2 /2m, a is the spin-orbit coupling con- 
stant, M = gfisMo-, and at (i = x, y, z) are the Pauli ma- 
trices. This Hamiltonian can be diagonalized easily, re- 
sulting in two eigen wave- functions <Ppp(r) = u A1 (p)e lp r 



and eigenvalues 



£/1 (p) =p72m +(-!)%, 



(2) 



(3) 



with X p = 



a 2 p 2 and p, = 1, 2. It is useful to in- 
troduce a unitary transformation U p = (ui(p), i*2(p)), 
by which the basis of the system is changed from a 
spin one to a helicity one. At the same time, Hamil- 
tonian becomes a diagonal matrix Hq = UpHoU p = 
diag(£i(p),e 2 (p))- 

We consider an electric current flowing along the x axis 
when the system is driven by a weak dc field E along 
the y-direction. In the helicity basis, the single-particle 
operator of this current, j x (p) = Up jx(p)U p , is given by 



3x{p) = e 



^(Xp-ma 2 ) ^(Mp x + iXp Py ) 
^L-{Mp x -iXp Pv ) -^-(Xp + ma 2 ) 



(4) 

and the corresponding macroscopical quantity is ob- 
tained by taking the statistical average over it, J x = 
S P Tr[j a; (p) j o(p)], with p(p) being the distribution func- 
tion. By definition, the anomalous Hall conductivity is 
determined by a xy = J x /E. 

In the spin basis, the lesser Green's function, G < , is a 
2x2 matrix and obeys the Dyson equation 



where U is the one-body external potential due to the 
dc field. The self-energies S < ' a arise from the electron- 
impurity interaction. In this paper, we consider a 
Coulomb interaction between electrons and impurities, 
which can be described by a potential V(p — k). This 
potential corresponds to a scattering of an electron from 
momentum state p to state k. Under homogeneous con- 
dition, Eq. (jSJ) in the momentum space can be written 
as 



■^ + <E-Y p \,>[p.T) -!-/[//,,/.(!>. T)} ■■■■ ~ ? 



dp 
dT 



scatt 

(6) 



with the Wigner distribution function p(p, T) = 
-iG<(p,T,T). 

We define a distribution function in helicity basis as 
p(p,T) — Up p(p,T)U p . To derive the kinetic equation 
for p, we multiply the Eq. © from left by Up and from 
right by Up. Due to the unitarity of the transformation, 
the right-hand side (rhs) of equation is obtained simply 
by replacing the Green's function in spin basis with that 
one in helicity basis 



dp_ 

dT 



dtil{X>,G<} + -{Z<,G>} + }(T,t>)(t\T). 

(7) 



At the same time, using the facts UpV p pU p — V p p — 

Vp[/+£/ p p-pl/+Vp[/p and V P U+U P = -U+V P U P , we 
find that the left-hand side (lhs) of the kinetic equation 
for distribution function in the helicity basis takes the 
form 



d 



dT 



E • V p \p - eE • [p, U+VpUp] + i[H , p}. (8) 



Further, we assume that the applied dc field is weak 
enough and only the stationary linear response of the 
system needs to be considered. Hence, we can divide 
the distribution function into two parts, p = po + 
pi, with unperturbed distribution function po(p) — 
diag(rip[£i(p)], np[e 2 (p)]) and the function in the first or- 
der of the dc field p\. 

To simplify relaxation term (jTJ, it is necessary to 
express the two-time lesser and greater Green's func- 
tions through one-time Wigner distribution. The pow- 
erful tool to perform this is the two-band generalized 
Kadanoff-Baym ansatz (GKBA), 22 which has been suc- 
cessfully applied to investigate the optical properties of 
semiconductors^ In the first order of dc- field strength, 
the GKBA reads 



Gf — — Ggpi + pi Go — GJpo + PoG", 



(9) 



[G 



-l 
o 



U, G < ] = XTG < + E < G a - G r S < — G < E a , (5) 



G> = -G r oPl + Pl G a + G{(I - po) - (/ - Po)G a 1 , (10) 

where, Gg° and G^' a are the Green's functions in the 
zero- and first-order of dc field, respectively. Note that 
G^' a are off-diagonal matrices. This fact is obvious 
from the Dyson equations which they should satisfy. In 
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Eqs. © and IjlOp. for shortness, the time arguments in 
Green's functions (iijia) are not written out. Under the 
stationary condition, the p\ is independent of the time. 
After these simplifications, the scattering term can be 

formally denoted as a sum of two components: ^ 



dT 



I 




II 


dp 


^ 


op 


scatt 


dT 


scatt 


dT 



II 



is off-diagonal and contains 

all the elements, associated with the functions G^' a . In 
this component, the perturbed distribution does not ap- 
pear because we are only interested in the linear response 



of the system. 



or 



includes all the remaining part. 



It is obvious that the driving force in the lhs of kinetic 
equation, i.e. the term proportional to E, comprises two 
components: the first of which, eE • V p poi contains only 
diagonal elements, while the another one is off-diagonal. 
According to this fact, we break the kinetic equation into 
two, 



eE- V p po + i[H ,p[} 



dp 

dT 



(11) 



distribution functions p\ and p{ 2 do not enter in the scat- 
tering term of off-diagonal Eq. I|12f> . Hence, this equation 

can be solved independently. On the other hand, the rhs 
1 



ofEq. (HU, $ 



generally depends on the distribu- 



tion p{ . However, as can be seen below, the contribution 



to Sfi- 
uo dT 



from the off-diagonal p{ is of higher-order of 

impurity density and can be neglected. Consequently, 
the Eqs. i(TT|) and (jT2|) become independent and can be 
resolved separately. 

We should note that there is a nonvanishing solution 
of the Eq. l|12f) if the collision term on the rhs is ignored. 
This pY leads to a Hall conductivity independent of any 
electron-impurity collision and results in the well-known 
intrinsic AHE. We will show below that neglect of the 
rhs of the Eq. I]12H is reasonable in the transport study. 
However, to resolve Eq. Qlljl . the scattering term can not 
be disregarded. Hence, the anomalous Hall conductiv- 
ity produced by the corresponding solution p\ becomes 
collision-related. 

B. Intrinsic anomalous Hall effect 



- eE • [p , U+V P U P ] + i[H , P ['] = - ^ 



(12) 



which have two solutions p\ and pi 1 , respectively: P \ 



In this paper, we restrict ourselves to consider the AHE 
being linear in dc field and to the leading order in the 
impurity-density expansion. In this case, the Eqs. 
and i|12|) are independent of each other. In fact, from 



the definition of 



11 



it is evident that the perturbed 



To carry out the expression of p^ 1 , we start with ana- 
lyzing the retarded and advanced Green's functions to 
the first order of dc field, G^' a , which appear in the 
scattering term of the rhs of Eq. (|12|l . It is well known 
that G^'° vanish in a one-band electron gasi^4 From the 
Dyson equations we can see that in the case of two 
bands, the diagonal elements of G^' a still vanish but their 
off-diagonal elements are nonzero and proportional to 
(GS' Q )n " (GS'°)22, G?° = ^eE • a z [G r ' a ,U+V p U p }. 
Inserting G^' a into the scattering term, we obtain 



(pfW(p) 



V P u+(p)up(p) / dujn F {uj) 



, M(G5W-(G5)a 

ep(p) -e n (p) 



(13) 



with p, = 3 — p. This expression can be further simpli- 
fied if we consider that the collision broadening of the 
perturbed Green's functions G ' a play secondary roles in 
transport and the imaginary parts of these functions re- 
duce to (5-functions. In result, the pi 1 has a simple form 



(14) 

Note that this result can also be derived if we simply ne- 
glect the scattering term on the rhs of Eq. 112(1. Taking 
the statistical average over the current operator j x (p), 
we find the contribution of p{ 7 to anomalous Hall con- 



ductivity: 



J 11 

■it) 



Ma 2 e 2 



dp 1 



(n F [ei(p)] - n F [e 2 (p)}). (15) 



This result agrees with the previous study on intrinsic 
anomalous Hall effect m& Obviously, it is independent of 
any electron-impurity scattering and connects with all 
the electron states under the Fermi surface. Further, it 
is clear that this expression of cr^ is related to the Berry 
phase and the topology of the energy bandsi^li Here, we 
will propose another interpretation of this intrinsic AHE 
by extending the treatment of Zhang and Yang in the 
study on spin-Hall effect^ 

In fact, the intrinsic contribution to AHE originates 
from a quantum interference between two bands per- 
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turbed by the external field and is associated with the 
nonvanishing interband dipole moment. The wavefunc- 
tion up to the first order of the electric field can be writ- 
ten as, 



\<P% > 



(16) 



probability to transit from one band to another. In result, 
the nonvanishing interband polarization emerges and re- 
lates to all the electron states. Note that there is no effect 
of this transition on diagonal elements of distribution up 
to the first order of dc field. 



with 



<^|eET|^> (0) ^ 

2^ ; tt~, l¥ Mk > 



■-tip 



Sak 



(17) 



being the first order perturbation. In order to take into 
account the statistical characters of the electrons, we will 
utilize the second quantization formalism. The field op- 
erators have the forms, ?/Vp = IVw> > c mp + IVw > c mp 
and ippfi — < ¥ , Mp| c ip+ < ^Mpl c jSp' with the electron cre- 
ation and annihilation operators, c^ p and c^ p . By def- 
inition, the interband polarization p^(p) describes the 
quantum interference of perturbed electrons in different 



bands, p MA (p) 
field, it reads 



< V'ppV'p/i >^ In the first order of dc 



(pi 1 )^(p) =< A« > { n F [e p (p)} - n K [ e „(p)]}, 



where the relations < 
and < c mp c w > 



(18) 

> 



^« >= - < 
n F[£p(p)] are used. Substituting Eq 



(|17|l into this equation and considering the fact that < 



m+ k eE 



i5(p — k), we finally 



arrive at expression i fTij l. 

It is obvious that this dc-field-induced transition can be 
interpreted as a linear stationary Rabi process. The well- 
known Rabi oscillation occurs in the presence of ac field 
and has been widely investigated from the viewpoints of 
optics, semiconductors, and atomic physics etc. However, 
the ac field is not the essential factor of occurring the 
interference, whereas it ensures the time oscillation of the 
Rabi process. If the applied ac field is replaced by a dc 
one, the quantum interference also may take place. For 
that, the necessary condition is the nonvanishing of Rabi 
frequency, which is proportional to the dipole moment, 
and the strength of external field. 

At the same time, from a perturbative point of view, 
this component of interband polarization can also be un- 
derstood as an interband transition between two unper- 
turbed bands. In this picture, all electrons have finite 



C. Disorder- mediated anomalous Hall effect 

The distribution functions p\ satisfy the Eq. . The 
lhs of this equation consists of a diagonal driving term 
eE • V p riF(e f j i (p)), and an off-diagonal matrix i[Ho,p{] 
depending only on off-diagonal elements of distribution 
function. Before seeking the solution of this equation, 
we should carry out the self-energy in the self-consistent 
Born approximation. It takes a complicated form and is 
presented in the Appendix. Due to the spin-orbit cou- 
pling, the wave functions possess an additional momen- 
tum dependence, leading to a nonglobal transformation. 
In result, each element of the self-energy, and hence the 
scattering term in Eq. (f 1 II) , becomes a function of all the 
elements of matrix Green's function. This fact indicates 
that the scattering can induce an admixture of the spin- 
orbit-coupled bands as well as the external field. In re- 
sult, some novel processes contributing to the distribu- 
tion functions appear. 

To simplify the relaxation term standing on the rhs 
of Eq. iJTTJ, we first analyze the lowest exponent in the 
impurity-density expansion for the elements of distribu- 
tion function. Note that the scattering always provides 
a contribution of order of iVj. Since the diagonal driv- 
ing term in the lhs of Eq. I|ll|) is independent of impu- 
rity density, the diagonal elements of distribution (p{) w 
should be of order of (Aj) -1 . We substitute these (p{) w 
into the off-diagonal parts of the relaxation term. Since 
the term i[Ho,p{] in the lhs of Eq. 111|) is proportional 
to the off-diagonal elements of the distribution, the lead- 
ing order of the off-diagonal (Px)ws in impurity-density 
expansion should be (A;) . Hence, all terms containing 



dp 

8T 



the off-diagonal (p{) M ^, as well as (p[ 

make contributions of higher-order of impurity density 
and therefore can be ignored. Under these considera- 
tions, we find the coupled equations for diagonal elements 
of distribution function, 



J 



-eE-Vp^FM^^I^p-k)! 2 ^!^ (19) 

k 

I 



with A M „ = £(£,u(p) -e v (k)), oi(p,k) 
a 2 kpcos(4> p - (f>k))/X p Xk and a 2 (p,k) 



(A p A fc + M 2 
(X p X k - M 2 



a 2 kpcos((j> p - <pk))/XpXk- 

The dc field is assumed to be applied along the y 
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axis. Hence, the lhs of Eq. l|19|l depends on the angle 
of momentum through a sine function. From the funda- 
mental triangle relation sin 4>k — ~ sin(0 p — </>k) cos <p v + 
cos(0 p — 0k) sin0 p and the symmetry argument that the 
terms with sin(0 p — </>k) vanish, it can be followed that 



the diagonal distribution has a simple angle-dependence 
[p[)vn = e E®fj.(p) sin0 p . Due to the elastic nature of 
electron-impurity scattering, functions can be car- 

ried out analytically 



dn F (£/x(p)) 



'2flp 



1 \ ggjjjg) 



dp 



-1 depjpn) 



liip, 



F 2mm) 



_-l ,-1 : 



(20) 



where p M is determined from equation, £p(p M ) = £ M (p), 
and the different relaxation times Ti^ v (i = l,.3,fi,v — 
1, 2) are defined by 

— = 27rn^|^(p-k)| 2 A VI ,(p,k), (21) 

T% ^ v k 

with Ai jUI/ (p,k) = i(l - cos(0 p - k ))ai(p,k)A MI/ , 
A 2 ^(p,k) = |a 2 (p,k)A M1/ and A 3jLU/ (p-k) = ±cos(</> p - 
0k)a 2 (p,k)A /iJ/ . 

From Eq. H19|) . we find that the terms in the rhs cor- 
respond two distinct scattering processes: the first two 
terms, i.e. the terms associated with A MM , describe the 
well-known intraband process, while the remaining come 



from an interband transition. It should be noted that 
the latter interband process entirely arises from the spin- 
orbit coupling: it results in a collision-related admixture 
of two bands. In the absence of spin-orbit interaction, as 
shown in the previous studies, p^p, at any time, does not 
appear in the equations for even when the external 
field is strong^ 

At the same time, we should note that ti mm is the longi- 
tudinal relaxation time defined in the framework of Bloch 
equations. However, the t 2ai „ and are novel relax- 
ation times and, up to now, have not been studied in the 
literature. 

Further, the off-diagonal elements of distribution are 
given by 



(p{)i2 = { P [)ti = ^ E i y (p - k )i 2ft 3(p> k)(-i)^{A w [(p[) MM ( P ) - (piwk)] - a^[(p{) m » - { P [) m m 

p k^=l,2 

with a 3 (p, k) = [akXp sin(0 p — </>k) + iaM(p — k cos(0 p — </>k))]/Afc A p . It is obvious that this function depends on the 
momentum angle not only through a sine, but also through a cosine function 

G°i)i2(p) = *Ci(p)sin0p + C2(p)cos0 p . (23) 

dip) (* = li 2) are determined by 
Ci(p) = TT ^ |^(p-k)|V(p,k)(-l)^{A w <f>4p)[l-cos(0p-0 k )]-A^[$4p)-a> p (p M )cos(0 p - ( /)k)]}, (24) 



C2(P)= 4A~ ^ inP-k)| 2 a 5 (p,k)(-l)nA w <I> M (p)sin(0p-0k)-A w - I <i> p (^)sin(0p-0k)} (25) 

p kp=l,2 

I 



with a 4 (p,k) = aM(p — fccos(0 p — 0k))]/A/cA p and 
as(p, k) = a 3 (p, k) - ia 4 (p, k). 

Similarly, from Eq. (1221) we can see there are two scat- 
tering processes to form the nonvanishing interband po- 
larization. One of which, associated with the terms pro- 
portional to A^£ , is the conventional interband transition 



of the perturbed electrons. Due to the scattering admix- 
ture of spin-orbit-coupled bands, even an intraband tran- 
sition, connecting with the first two terms of expressions 
(|24|l and l|25|l , makes a contribution to the interband po- 
larization. 

Based on these arguments about the scattering pro- 
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cesses in Eas. I|19[l and ()22fl. we can generalize our re- 
sults to the case of arbitrary two-band 2D system. In 
a two-band system with noninteracting wave function 
u fl (p)c ip ' r , the intcrband or intraband processes, in- 
duced by impurity scattering, are described through 



u+(p)V(p — k)u„(k). If the off-diagonal distribution 
functions in the scattering term can be omitted, by means 

of the well-known Fermi golden rules, the ^ can be 
written as 



dp_ 

dT 



scatt, /j, v 



= nJ2\V( P - k)| 2 U +(pK(k)u+(kK(p) 



x {A^[(pf) w ( P ) - (piw(k)] + am^p) - (Aum} 

I 



(26) 



Note that this equation holds for a general two-band 2D 
system. 

After distribution functions are obtained and the sta- 
tistical average over j x is taken, the anomalous Hall cur- 
rent, related to the longitudinal transport, can be deter- 
mined by 



4 = 2£ 



uMt 2 {v) 

Xr, 



■ cos 



a(i (p) sin 2 



(27) 



Since the diagonal elements of the current operator is 
proportional to p x , the contribution of diagonal distribu- 
tions to anomalous Hall current vanishes. It is obvious 
that such anomalous Hall current has a magnitude of or- 
der of (TV,;) . 

The procedure for deriving the off-diagonal distri- 
bution tells us that this collision-related contribution 
to anomalous Hall effect originates from a disorder- 
mediated process. The electrons, influenced by exter- 
nal dc field, experience the impurity collision and join in 
longitudinal transport. Such scattering produces diago- 
nal elements of distribution function of order of (JVj) -1 . 
At the same time, these electrons scattered by impurities 
give rise to an interband polarization of order of (iVj)°. 
The disorder only plays an intermediate role. Note that 
at each stage of this disorder-mediated process, both in- 
terband and intraband scatterings occur. 

Note that this disorder-mediated mechanism and the 
well-studied side-jump mechanism^ are identical physi- 
cally. Contributions to Hall conductivity from both the 
mechanisms are independent of the impurity density but 
rely on the electron-disorder collision. At the same time, 
they all are related to the electron states near the Fermi 
surface. However, formally, these two mechanisms are 
completely different. It is well known that the side-jump 
process corresponds to a lateral displacement of the cen- 
ter of the wave-packet during the scattering when the 
spin-orbit interaction is included into the potential of the 
electron-impurity scattering. It is associated with the 
scattering-related term in the current operator*! How- 
ever, in our study for Rashba spin-orbit coupling, the 
current operator becomes impurity-independent. 



III. NUMERICAL RESULTS 

We also perform a numerical calculation for a 
GaAs/AlGaAs heterojunction to investigate the intrinsic 
and disorder-mediated anomalous Hall effect in a Rashba 
two-dimensional electron system with magnetization. In 
calculation, we first consider a long-range scattering be- 
tween electrons and the remote impurities separated at a 
distance s. The corresponding potential takes the form: 
V(p) = U(p)/ K (p) with U( P ) = ^e-*I(p)» Hp) = 
fo 3 (b + p)~ 3 is the form factor. The wave-function param- 
eter b is given by 6 3 = f^#(A + ff A dcp ). N c and A dcp 
are the electron density and the density of depletion layer 
charges, respectively. n(p) = 1 + q s H(p)/p is the factor 

2 

coming from the Coulomb screening. q s = "' e and 
H{p) = i(86 3 + 9pb 2 + 3bq 2 )(b + q)~ 3 . In calculation, we 
take m = 0.067m o , s = 500 A and A dcp = 1 x 10 13 m~ 2 . 
In our formalism, the intrinsic anomalous Hall effect can 
be obtained directly by evaluating integral (|15|l . To carry 
out the disorder-mediated anomalous Hall current by Eq. 
(I27|l . we should determine first the diagonal elements of 
the distribution from Eqs. H19fl and 120|) and then the off- 
diagonal ones by Eq. (|23[l . For long-range disorders, the 
results are plotted in Fig. 1 and 2. 

In Fig. 1, we plot the disorder-mediated a xy , intrinsic 
a xy and total anomalous Hall conductivity <j xy as func- 
tions of the spin-orbit coupling constants. It can be seen 
that the values of a xy and a~. are comparable. For dif- 
ferent magnetization M, with increasing the SO coupling 
constant, the intrinsic <r*.y always increases, while the a xy 
exhibits a complicated behavior: for large M it always 
increases as well as a xy , but for small M it first increases 
and then falls, and even becomes negative for large a. 
The similar behaviors can also be seen for total a xy . It is 
clear that a xy and a xy do not always have the same sign 
and a compensation occurs for small M and large a. 

We also calculate the dependence of anomalous Hall 
effect on the magnetization for different SO coupling con- 
stants. The results are shown in Fig. 2. Again, we see a 
sign change for small M and large a. Besides, with in- 
creasing the magnetization above a critical value, the a xy 
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FIG. 1: Dependencies of disorder-mediated a xy , intrinsic ai y 
and total anomalous Hall conductivity a xy on the spin-orbit 
coupling constant in a Rashba 2D electron system with mag- 
netization. The lattice temperature T = K and the electron 
density 7V C = 1 x 10 11 cm" 2 . 

falls. Taking one with another, in the studied parameter 
range, the absolute values of disorder-mediated a^ y are 
always lesser than the intrinsic a xy . 

Further, to demonstrate the collision-related feature 
of the a xy , we also compute the Hall conductivity for 
short-range collision by replacing the potential V(p — k) 
by a momentum-independent one u. Note that the re- 
sultant a^ y is independent of the magnitude of u. The 
calculated results are plotted in Fig. 3. It can be seen 
that the a T xy from the short- and long-range collisions ex- 
hibit completely different behaviors in the dependencies 
of the SO coupling constant and magnetization. In the 
parameter regime: a < 2 x 10 _9 eVcm and M < 3meV, 
where the long-range disorders have strong effect on a xyi 
the Hall conductivity practically can not be affected by 
the short-range disorder. At the same time, in Fig. 3(b) 
an abrupt decrease of a xy with increasing the M is seen. 
Both these features can be understood from the fact that 



12 3 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
M (meV) 



FIG. 2: Disorder-mediated ai y , intrinsic a xy and total 
anomalous Hall conductivity a xy as functions of the magne- 
tization. The other parameters are the same as that in the 
Fig. 1. 



for short-range collision, the disorder-mediated mecha- 
nism has effect on Hall conductivity approximately when 
M becomes larger than the chemical potential. Note that 
this fact agrees with the study in Ref.l^iJ In Fig. 3(a), 
the a xy gradually descends with ascending a for large M 
since these M are always larger than the chemical poten- 
tial in the studied regime of a. In contrast to the case of 
long-range scattering, disorder-mediated contribution in- 
duced by the short-range collision always has an opposite 
sign with respect to the intrinsic one. 

We should note that in the discussion above, to exhibit 
the effect of different mechanisms on AHE, we divide the 
contribution to Hall conductivity into an intrinsic one 
and a disorder-mediated one. Actually, in experiment, 
both can not be distinguished. They together result in a 
measurable Hall resistivity proportional to the square of 
the impurity density. 
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a(10" 9 eV cm) 
12 3 4 




M (meV) 

FIG. 3: Different effects of short- and long-range disorders 
on the anomalous Hall conductivity. The coupling-constant 
and magnetization dependencies of total a xy for short-range 
collision are plotted respectively in the insets of Fig. 3(a) and 
Fig. 3(b). 



IV. CONCLUSION 

We have proposed a two-band kinetic equation ap- 
proach to the anomalous Hall effect in a 2D elec- 
tron system with a Rashba spin-orbit interaction and 
an exchange-field-induced magnetization. The obtained 
equation has been resolved by considering the electron- 
impurity collision in the self-consistent Born approxi- 
mation, ft is clear that there exist two contributions 
to the impurity-density-free AHE. One of which arises 
from a dc-field-induced transition between two spin- 
orbit-coupled bands, which can also be understood as 
a linear stationary Rabi process. Another contribution 
is related to the collision, but is independent of the im- 
purity density. It comes from a process mediated by dis- 
order: electrons participating in longitudinal transport 
are scattered again by impurity, yielding a nonvanish- 
ing intcrband polarization. Numerically, we have demon- 
strated the dependencies of both these contributions on 
the magnetization and SO coupling-constant : they al- 
ways exhibit a compensation for short-range disorder, 
while such compensation occurs in the case of long-range 
collision only for small magnetization and large SO cou- 
pling constants. 

Since obtained Eqs. (|14|l and (|26[1 are valid even for a 



general two-band 2D system, the interband polarization 
consisting of such two terms is a universal property of 
2D systems with spin-orbit interaction. In result, the 
quantities, connecting with the interband polarization, 
such as anomalous Hall current, spin-Hall current etc. 
should be formed from two distinct processes, which are 
associated with the electron states under and near the 
Fermi surface, respectively. 
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APPENDIX: SELF-ENERGIES 

In the spin basis, electrons experience a spin- 
independent long-range disorder collision, which can be 
described by a potential V^(p). The self-energies for re- 
tarded, advanced and lesser Green's functions in the self- 
consistent Born approximation take the forms 

£^<(p) = J2 \V(p - k)| 2 G^<(k). (A.l) 

k 

When a transformation from spin basis to hclicity ba- 
sis is performed, the forms of the self-energies £ r > a >< = 
U+W' a '<U p are changed to 

E r >°><(p) = |V(p-k| 2 C/+[/ k G^ < (k)f/+[/ p . (A.2) 

k 

It can be seen that, due to the non-global feature of 
the transformation, the scalar scattering potential in 
spin basis becomes a matrix in hclicity basis T(p,k) = 
Up V(p — k)[/k- Generally, each element of self-energy 
matrices contains all the elements of Green's functions. 

It is convenient to express the self-energies through a 
sum of several matrices. Therefor, we first define several 
matrices: diagonal matrices B^ (i = 1..4) 

(Bi) w = ^(G + a z Ga z ) w = G w , (A.3) 



(Ba)^ = -{a x Ga x + (JyGcjy)^^ — Gfifi, (A. 4) 



(Ba) w = ^([o'x, G] + <r z [a xl G]a z )^ = G AAI - G^, 

(A.5) 
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(B/Owi = 2 ([°V) G]+o\j[o- y , G]cr z ) w = (-l) A 'i(G Al p+Gp A1 ), (Ca)^ = ^(E *) ^] ~ ^I ^ G]a z )^ = Gpp - Gp M , 

(A.6) " (A.9) 

off-diagonal matrices Gj (i = 1..3) and matrices = o^C*, (Dj)^ = (— l) M+1 (Ci) A1 p. By 

means of these matrices, we can rewrite the self-energies 
as 



(Ci) M p = -(G - a z Qu z )^ = G MA , (A. 7) 



(02)^/2 — ^(^Gcx — PyGdy)^ — Gfifj,, (A. 



S(p)= l^(p-k)| 2 6,:B l (k)+ I^P-kjriciCiCkJ + diDiCk)]. (A.10) 

k,i=1..4 k,z=1..3 

&i (i = 1..4), c, and (i = 1..3) are the factors depending on both the momenta p and k, 

' -[M 2 -(-l) l A p A fc + a 2 fcpcos(0 p -0 k )],i = l,2 (A.ll) 



h = — sin(>p - k ), 64 = -— — [k -pcos(0 p - k )], (A.12) 

ZAp ZApAk 



2 kp+{M 2 + {-l) j X p X k ) cos(<£ p - <£ k )],j = 1,2 (A.13) 



%otk iJVl 1 

c 3 = — sin(0 p - k ), = ^^(A P - (-l) J A fc ), j — 1,2 d 3 = ——(akMcos(<f> p - k ) - apM). (A.14) 

ZAfc ZAkAp ZApAk 

In comparison to the systems without magnetizationj2I the self-energies take more complicated forms due to the 
additional dependence of wave functions on the module of momentum. 
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